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We have used the COBE DIRBE database to derive near- and mid-infrared light curves for a 
well-defined sample of 38 infrared-bright Mira variable stars, and compared with optical data from 
the AAVSO. In general, the 3.5 /im and 4.9 /im DIRBE bandpasses provide the best S/N light 
curves, with S/N decreasing with wavelength at longer wavelengths. At 25 /im, good light curves 
are only available for ~10% of our stars, and at wavelengths > 60 /im, extracting high quality 
light curves is not possible. The amplitude of variability is typically less in the near-infrared than 
in the optical, and less in the mid-infrared than in the near-infrared, with decreasing amplitude 
with increasing wavelength. On average, there are 0.20 ± 0.01 magnitudes variation at 1.25 /im 
and 0.14 ± 0.01 magnitudes variation at 4.9 /im for each magnitude variation in V. The observed 
amplitudes are consistent with results of recent theoretical models of circumstellar dust shells 
around Mira variables. For a few stars in our sample, we find clear evidence of time lags between 
the optical and near-infrared maxima of phase ~ 0.05 — 0.13, with no lags in the minima. For 
three stars, mid-infrared maximum appears to occur slightly before that in the near-infrared, but 
after optical maximum. We find three examples of secondary maxima in the rising portions of 
the DIRBE light curves, all of which have optical counterparts in the AAVSO data, supporting 
the hypothesis that they are due to shocks rather than newly-formed dust layers. We find no 
conclusive evidence for rapid (hours to days) variations in the infrared brightnesses of these stars. 

Subject headings: stars: AGB and post-AGB, stars: variable: Miras 
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1. Introduction 

Current theories of stellar evolution indicate 
that low mass (<8 M Q ) stars pass through the pul- 
sating Mira star phase at the end of their time as 
Asymptotic Giant Branch stars (Wood 1997; Olof- 
sson 1999). Although the time that stars spend as 
Mira variables is relatively short (^6 x 10 4 years; 
Wood (1990)), this stage is very important be- 
cause these stars have very high rates of mass loss 
(up to 1CT 4 Mq yr" 1 ; Knapp & Morris (1985)). 
This means that they return large quantities of 
carbon, oxygen, and nitrogen-enriched gas to the 
interstellar medium (Sedlmayr 1994). At present, 
the mass loss process, the properties of the cir- 
cumstellar envelopes of these stars, and how these 
relate to the variability of the star are still not very 
well-understood. 

To better understand these processes, detailed 
studies of Mira variables at infrared wavelengths 
are helpful. In the near-infrared, photosphcric 
emission dominates, while at wavelengths longer 
than ~3 /im emission from dust in the circum- 
stcllar envelope becomes important (Woolf & Ney 
1969). Infrared observations have been used to 
constrain models of the circumstellar envelopes 
of Mira variables and to estimate mass loss rates 
from these stars (Jura 1986; Bedijn 1987; Onaka 
ct al. 1989a, b; Anandarao et al. 1993; Le Bertre & 
Winters 1998). These studies, however, are based 
on time-averaged broadband fluxes from the In- 
frared Astronomical Satellite (IRAS) (Jura 1986; 
Bedijn 1987; Anandarao et al. 1993; Le Bertre & 
Winters 1998) or time-averaged 8 — 23 (im IRAS 
Low Resolution Spectrometer (LRS) spectra (On- 
aka et al. 1989a,b), and have not investigated mid- 
or far-infrared variability of these stars. 

To better understand mass loss in these stars, 
infrared monitoring is helpful. Such studies have 
revealed variations in the infrared spectral en- 
ergy distributions of Miras with time. Broad- 
band infrared observations indicate that the am- 
plitude of variation tends to decrease with wave- 
length from the optical through the mid-infrared 
(Lockwood & Wing 1971; Harvey et al. 1974; 
Catchpolc et al. 1979; Barthes, Chencvez, & Mat- 
tel 1996; Le Bertre 1992, 1993; Little-Marenin ct 
al. 1996). Spectroscopic observations with IRAS 
(Little-Marenin & Stencel 1992; Hron & Aringcr 
1994; Little-Marenin ct al. 1996), ground-based 



telescopes (Creech-Eakman ct al. 1997; Monnicr 
ct al. 1998), and the Infrared Space Observatory 
(ISO) (Onaka et al. 1997) showed that at least 
for some Miras the silicate features tend to de- 
crease in flux from light curve maximum to min- 
imum. There is also evidence from IRAS data 
for far-infrared variations of Miras (Little-Marenin 
& Stencel 1992; Little-Marenin et al. 1996; Spear 
ct al. 1993), however, the IRAS light curves are 
poorly-sampled, so variations at these wavelengths 
are not well-studied. 

Infrared monitoring of Mira variables has also 
revealed wavelength-dependent phase lags. The 
maxima at 1 /im (Petit & Nicholson 1933; Lock- 
wood & Wing 1971; Barnes 1973) and 2.7 
(Maran et al. 1977) typically lag those at V by 
about 0.1 phase, while no lags are seen in the min- 
ima. At present, phase lags at longer wavelengths 
are not well-determined. There have also been a 
few reported observations of large- amplitude (0.5 
— 1 magnitudes) very short-term (hours to days) 
variations in the visual (Maffci & Tosti 1995; de 
Laverny et al. 1998) and near-infrared (Smak & 
Wing 1979b; Guenter & Henson 2001) light curves 
of Mira variables, but such variations have not yet 
been reported at longer wavelengths. 

The previously-cited studies generally provide 
long-wavelength infrared light curves for only a 
relatively small sample of stars, and generally sam- 
ple only a fraction of the light curve of these stars. 
To obtain more complete infrared light curves for 
a larger sample of Mira stars, we have used data 
from the Diffuse Infrared Background Experiment 
(DIRBE) instrument (Hauser et al. 1998) on the 
Cosmic Background Explorer (COBE) spacecraft 
(Boggess et al. 1992). To date, little research has 
been done on infrared variability of stars using 
DIRBE data, in spite of its suitability for such 
studies. 

The COBE satellite was launched in Novem- 
ber 1989, and the DIRBE instrument functioned 
at cryogenic temperatures from 11 December 1989 
until 21 September 1990. This 10 month timescale 
is suitable for searching for variations in long pe- 
riod variables. In the course of a week, a typ- 
ical point on the sky was observed 10-15 times 
by DIRBE, and over the course of the cryogenic 
mission, it was observed about 200 times (Hauser 
ct al. 1998). The temporal coverage of DIRBE 
varies with sky position. Some parts of the sky 
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had relatively complete coverage; for other re- 
gions, there were gaps in the coverage of several 
months (Hauser et al. 1998). DIRBE operated at 
10 infrared wavelengths (1.25, 2.2, 3.5, 4.9, 12, 
25, 60, 100, 140, and 240 /im), providing a nearly 
complete infrared spectrum. This allows us to tie 
stellar phenomena, as seen in the visual and near- 
infrared, to circumstellar emission, seen at longer 
wavelengths. The beam size of the DIRBE in- 
strument, 0.7 degrees, although large, is still suffi- 
ciently small that individual infrared-bright stars 
can be discerned. 

We have used DIRBE data to extract 1.25 nm 
— 240 fim light curves for 38 Mira variables and 
compared them to available optical data. We have 
investigated the shapes of the light curves, ampli- 
tude variations, and phase lags, and have com- 
pared these results to recent models of circumstel- 
lar dust emission from Miras. 

2. Sample Selection and Infrared Light 
Curve Derivation 

Our stars were selected from the 12 ^im flux- 
limited sample of variable stars described by Sloan 
& Price (1998) and Sloan et al. (1998). This sam- 
ple was obtained by cross-referencing the General 
Catalogue of Variable Stars (Kholopov et al. 1985- 
1988) with the IRAS Point Source Catalog (1988). 
We selected a subset of this sample, with an IRAS 
12 fxm flux density cutoff of 235 Jy and classifica- 
tion as Mira stars. A total of 38 stars are included 
in our set (Tables 1 and 2). Eight of these are 
classified as carbon-rich Mira variables (Table 1), 
while 30 are oxygen- rich (Table 2). In Tables 1 and 
2, we provide the positions, the period of variabil- 
ity, the IRAS 12 /im magnitude, the optical spec- 
tral type, and the dust spectral type (based on 
the IRAS LRS spectra). The positions are from 
the SIMBAD database; the rest of the informa- 
tion is from Sloan & Price (1998) or Sloan et al. 
(1998). 

Infrared light curves were derived for these stars 
from the DIRBE data using two methods. We 
first extracted light curves from the COBE DIRBE 
Weekly Sky Maps. These are weekly- averaged in- 
tensity maps of the sky in the ten DIRBE wave- 
length bands. The stellar photometry was deter- 
mined by multiplying the surface brightess of the 
Weekly Sky Maps at the stellar position by the 



DIRBE beam solid angle, after sky subtraction. 
For sky values, we used the average of the fluxes 
at four positions 2° N, S, E, and W in equatorial 
coordinates. Flux uncertainties were obtained by 
using the quoted DIRBE photometric uncertain- 
ties combined in quadrature with the dispersion 
in the four sky measurements. At 1.25 /im — 240 
fim, respectively, the average noise levels in these 
light curves are 35 Jy, 40 Jy, 35 Jy, 20 Jy, 120 Jy, 
155 Jy, 140 Jy, 460 Jy, 2850 Jy, and 1250 Jy. 

The second method of obtaining light curves 
used the DIRBE Calibrated Individual Observa- 
tions (CIO) data files. The CIO database con- 
tains the calibrated individual l/8th second sam- 
ples taken in science-survey mode during each day 
of the cryogenic mission. For all scans that passed 
within 0.3° of the target position, a linear base- 
line was fit to the sections ±1.35° — 2.25° from 
the point of closest approach. The point source 
photometry was obtained by subtraction of this 
baseline and correcting for DIRBE beam response. 
The uncertainties in the point source photometry 
are calculated as the quadrature sum of the rms 
noise of the baseline, an error due to positional 
uncertainties of V in both the in-scan and cross- 
scan directions, an error due to short-term detec- 
tor gain variations, and signal-dependent detector 
noise. For the light curves derived from the CIO 
data, average noise levels per datapoint are 25 Jy, 
20 Jy, 20 Jy, 10 Jy, 30 Jy, 55 Jy, 320 Jy, 765 Jy, 
4800 Jy, and 2750 Jy for the DIRBE 1.25 - 240 
/im filters, respectively. 

Datapoints with large uncertainties (greater 
than three times the average uncertainty) were 
eliminated from the dataset, since anomalously 
large errorbars were likely caused by cosmic ray 
hits or the presence of a nearby companion in the 
DIRBE scan. Careful comparison of individual 
DIRBE scan locations with the IRAS Point Source 
Catalog (1988) and the all-sky near-infrared cat- 
alog developed by Smith (1995) for position re- 
construction for the Infrared Telescope in Space 
(IRTS) mission (Murakami et al. 1996; Wheaton 
et al. 1997) shows that when an infrared-bright 
companion was present in the baseline portion 
of the scan, the derived stellar photometry was 
sometimes discrepant with large uncertainties. 

These two methods provided consistent light 
curves, with flux densities that agree within 15%, 
4%, 12%, 2%, 6%, 3%, and 12%, for the 1.25 
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/im - 60 /zm bands, respectively, when averaged 
over our light curves. As noted in the DIRBE Ex- 
planatory Supplement (Hauser et al. 1998), point 
source photometry from the Weekly Sky Maps 
may suffer a small systematic offset when com- 
pared to the CIO data. The CIO data are believed 
to be more accurate for point source photometry 
(Hauser ct al. 1998), so for the following analy- 
sis we focus on the CIO light curves. Note that 
for the six shortest wavelengths, the uncertainties 
in the CIO light curves are less than those of the 
photometry derived from the Weekly Sky Maps, in 
spite of the shorter time coverage per data point. 
This provides a second reason for using the CIO 
light curves rather than the light curves from the 
Weekly Sky Maps. 

Note that the uncertainties used in this pa- 
per do not include absolute flux calibration er- 
rors, and that color corrections were not made 
to the DIRBE flux densities. The color correc- 
tions are typically quite small, particularly at the 
shorter DIRBE wavelengths. For the stars in our 
sample, the color correction tables in Hauser et 
al. (1998) give estimated color correction factors 
(F 'observed /F true) that range from 0.99 - 1.03 at 
1.25 /zm, 0.99 - 1.03 at 2.2 /zm, 0.98 - 1.01 
at 3.5 /zm, 0.99 - 1.01 at 4.9 /zm, 0.89 - 1.1 
at 12 /im, and 0.5 — 2.8 at 25 fim. The range 
quoted for each wavelength reflects the uncertainty 
associated with the shape of the emissivity law 
(e oc v 1 ^ 2 ), the fact that the spectral energy dis- 
tribution of each star is variable in time, and the 
fact that each star has a slightly different spectral 
energy distribution. 

3. Optical Data 

For the period of the COBE mission, 23 out 
of our 38 stars have well-sampled visual (V) light 
curves available from the American Association of 
Variable Star Observers (AAVSO) dataset (Mattei 
2000). These stars are identified in the last column 
of Tables 1 and 2. These optical light curves are 
plotted along with the infrared data in Figures 1 
and 2. 



4. Results 

4.1. Overview of the Light Curves 

In Figure 1 we present plots of flux density vs. 
time for the six shortest wavelength DIRBE bands 
for six stars with relatively complete DIRBE light 
curves. In Figure 1, we also include the visual 
light curves when available, as well as the 2.2 /im 
to 12 /zm flux ratio as a function of time. For RW 
Vel, which does not have optical data available, we 
also plot F4.g /Jm /F2.2 p,m vs. time. In Figure 2, we 
provide 4.9 /zm light curves for the remaining 31 
stars in our sample, along with optical light curves 
when available. These same data are displayed 
in a different manner in Figure 3, where we plot 
the optical-infrared spectral energy distributions 
for all 38 stars at one or more times during the 
COBE mission vs. wavelength. 

In general, the 3.5 /zm and 4.9 /im light curves 
have the best S/N. At 3.5 /im and 4.9 /im, re- 
spectively, on average the stars have 18.0 and 17.1 
5d weekly-averaged flux measurements from the 
Weekly Sky Maps (out of a total possible of 41 
weeks). In the CIO light curves, at 3.5 /im, the 
average number of datapoints is 490, and at 4.9 
/im, 500 points are available on average. Almost 
complete 4.9 /im light curves (>30 weeks of 5a 
data) are available for five stars in the sample, 
while 31 stars (82% of the sample) have at least 
10 5tr datapoints available. Beyond 4.9 /im, the 
S/N drops with increasing wavelength. At 12 /zm, 
for example, out of the 38 sample stars only 14 
(37%) have more than ten datapoints detected at 
the 5(7 level or better in the Weekly Maps. By 25 
/zm, this number is down to nine stars (24%), and 
at the four longest DIRBE wavelengths, no stars 
have more than ten 5a weekly datapoints in their 
Weekly Map light curves. 

The stars in our sample are expected to be 
among the brightest in the infrared, so these statis- 
tics show the limitations of DIRBE data in study- 
ing infrared variability of unresolved sources. At 
60 /tm and longer wavelengths, the signal-to-noise 
ratio is too poor to obtain light curves using 
DIRBE data. At 25 /im, good light curves (20 
weeks at 5a) are available for only four out of our 
38 stars (11%). 
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4.2. Infrared Variability of our Sample 
Stars from DIRBE 

In many of the light curves shown in Figures 
1 and 2, variations are clearly apparent in the 
DIRBE data. Variations are also clearly evident 
in the spectral energy distributions plotted in Fig- 
ure 3. The best wavelengths at which to observe 
variability are 3.5 /mi and 4.9 /mi, where 30 and 
31 stars, respectively, show >3cr variability in the 
Weekly Maps. At 25 /um, only 3 stars vary more 
than 3cr in the weekly data. 

As an approximation based on the general ap- 
pearance of the light curves in Figure 1, when 
possible we fit the infrared light curves (in mag- 
nitude units) to sine functions to estimate ampli- 
tudes and dates of maxima (Table 3). Such fits 
were only possible for three stars in our sample: 
R Hor (Figure la), RW Vel (Figure lb), and R 
Cas (Figure lc). The rest of the light curves were 
too incomplete or too peculiar to be accurately fit 
to a sine function. Even for R Hor, RW Vel, and 
R Cas, these fits should be taken as approxima- 
tions, since the light curves are not perfect sine 
functions. These fits, however, provide useful es- 
timates of the amplitudes of variation and the JD 
dates of maxima as a function of wavelength (Ta- 
ble 3). Since we have less than a full pulsation cy- 
cle of infrared data for these stars, we were unable 
to obtain a good fit for the period of variability 
from the DIRBE data alone. We therefore fixed 
the period to be equal to that in the optical, and 
fit the DIRBE data for the amplitude and time 
of maximum. For R Hor and R Cas, we obtained 
the optical period by fitting the AAVSO data for 
JD = 2447500 - 2448700 (-3 pulsation cycles) to 
sine functions. We found periods of 416.5 days and 
439.7 days for R Hor and R Cas, respectively, close 
to the values given in Sloan & Price (1998) (Table 
2). For RW Vol, which does not have AAVSO data 
available for this time period, we simply used the 
period given in Sloan & Price (1998) of 443.1 days 
(Table 2). 

The stars in our sample with almost complete 
DIRBE light curves (Figure 1) can be divided 
into two groups: those with approximately sinu- 
soidal infrared light curves (R Hor, R Cas, and 
RW Vel), and those which show a plateau, or in- 
flection point, in the rising portion of their light 
curves (R Car, T Cep, and S Cep). Note that all 



three of the stars with sinusoidal light curves are 
oxygen-rich stars, while two stars with inflection 
points are oxygen-rich (R Car and T Cep) and one 
(S Cep) is a carbon star. 

The best example of a secondary maximum is 
seen in the light curves of R Car (Figure lb) at 
JD - 2448045 (phase - 0.75), where it is visible 
at all wavelengths between 1.25 [im and 12 /im, 
and maybe at 25 /im as well, with a duration of ~ 
40 days (Figure 1). These secondary maxima show 
an increase in brightness of ^15 — 40% relative to 
that expected from a sinusoidal light curve. This 
inflection point is also visible in the optical light 
curve, at a similar level (^0.4 magnitudes). 

The lower S/N light curves of T Cep (Figure 
lc) show weak inflection points at JD ~ 2448020 
at 1.25 /im, 2.2 /im, 3.5 /im, and 4.9 /im, with a du- 
ration of ~ 90 days. S Cep (Figure la) has possible 
inflection points between JD ~ 2447980-2448020 
in the four shortest wavelength DIRBE bands. In 
both cases, optical counterparts are indicated by 
the AAVSO data, with Amag(V) - 0.4. 

We also searched the CIO light curves for large 
amplitude, short-term (2 hours to 10 days) varia- 
tions in brightness, but found no conclusive de- 
tections greater than 0.2 magnitudes at the 3<r 
level. A handful of datapoints with very dis- 
crepant fluxes were noted in the CIO light curves, 
however, these were always individual isolated 
points, with discrepant fluxes at only one wave- 
length, thus their anomalous fluxes are likely due 
to cosmic ray hits rather than stellar variations. 
We note that in the CIO data, the typical time 
between successive observations of the same star 
is ^9 hours, so these data are not particularly sen- 
sitive to extremely short variations. 

4.3. Wavelength-Dependent Phase Lags 

We used our light curves to search for phase lags 
of the maxima as a function of wavelength. Five 
of our stars (R Hor, R Cas, R Cen, R Aql, and 
R Aqr), all of which are oxygen- rich, have both 
optical and infrared data available near maximum 
light. All five stars show hints of an infrared delay 
relative to V. 

The best example is R Hor (Figure la), which 
clearly shows a delay of the infrared relative to 
V at 1.25 /xm, 2.2 /im, 3.5 /im, and 4.9 /im, and 
perhaps at 12 /im and 25 /xm as well. The opti- 
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cal light curve peaks at JD ~ 2447930, while the 
infrared light curves are still rising until at least 
JD ~ 2447960, after which there are gaps in the 
DIRBE data until JD - 2447990. These gaps are 
followed by a decrease in the infrared brightness. 
This phase lag can also be seen in Figure 3d, the 
spectral energy distribution plots. The optical flux 
was brighter in the JD 2447897 spectrum than in 
the JD 2447950 spectrum, while the infrared fluxes 
were brighter at JD 2447950. 

The existence of this lag is supported by our 
sinusoidal fits to the light curves (Table 3), which 
show that the 1.25 /im maximum followed that 
at V by —53 days (0.13 phase) and the 2.2 /zm 
maximum lagged V by —65 days (0.16 phase). In 
Figure 4a, for the three stars for which we obtained 
sinusoidal fits to the light curves, we plot the phase 
lag relative to V as a function of wavelength. For 
RW Vel, for which we have no optical data, we 
assumed a 1.25 /im lag relative to V of 0.1 phase. 

R Cas also shows evidence for an infrared to 
V time lag (Figure lc; Figure 4). In the visible, 
R Cas was brightest at JD - 2447930, while the 
1.25 /xm — 25 /im emission was still increasing 10 
days later, at JD - 2447940 (phase - 0.02), after 
which there is a gap in the DIRBE coverage. The 
best-fit sine curves (Table 3) show a 1.25 /im lag 
relative to V of 38 days (0.09 phase) and a 2.2 /im 
lag relative to V of 47 days (0.11 phase). 

Our data suggest not just offsets between op- 
tical to infrared maxima, but also offsets between 
the DIRBE wavelengths. For two of our stars with 
relatively complete, high S/N DIRBE light curves, 
R Hor and R Cas, the shapes of the observed light 
curves differ from wavelength to wavelength. In 
particular, comparison of the 2.2 /im, 4.9 /im, and 
12 /im light curves hint that the 4.9 /im and 12 
/im maxima for these stars may have occurred be- 
fore the 2.2 /im peaks (but after the optical max- 
ima). The existence of these offsets is supported 
by our sinusoidal fits to the data (Table 3; Figure 
4). For R Hor, the best fit maxima at 1.25 /im to 
3.5 /im arc at JD = 2447993 ± 6 days (0.13 - 0.16 
phase), while those at 4.9 /im to 25 /im were 30 
— 50 days earlier (0.06 - 0.12 phase). For R Cas, 
the maxima at 1.25 /im to 3.5 /im occurred at JD 
= 2447999 ± 6 (-0.10 phase), while the 4.9 /im 
to 25 /im peak is at JD = 2447979 ± 3, 20 days 
earlier (—0.05 phase). 

Of the stars with no optical data available, the 



most complete infrared light curves are available 
for RW Vel (Figure lb). For this star, the max- 
ima at 1.25 /im, 2.2 /im, and 3.5 /im occurred at 
JD = 2448023 ± 9, while the 4.9 /im maximum 
happened 40 days earlier (0.09 phase) (Table 3). 
This is also visible in the F4.9 lim /^ > 2.2 time 
plot (Figure lb). Other stars with possible phase 
lags include R Cen, R Aqr, and R Aql, whose 1.25 
/im maxima appear to lag those at V by —30 days 
(0.05 — 0.1 phase). These are uncertain, however, 
because of incompleteness in the light curves. 

In addition to a mid-infrared/near-infrarcd 
phase shift, there may be smaller offsets between 
the three near-infrared bands. In particular, note 
that for the three stars with sinusoidal fits to their 
light curves (Table 3; Figure 4), the 2.2 /im lags 
are always larger than those at 1.25 /im and 3.5 
/im. This result, however, is uncertain because 
these light curves are not perfect sine curves. 

For six stars in our sample, there are both op- 
tical and infrared data available for a minimum in 
the light curve: R Cen, V Cyg, T Cep, R Lep, 
R Aqr, and R Car. In all of these cases, there is 
no evidence for a wavelength-dependent time de- 
lay. Of the six, the highest S/N data are those 
of R Car (Figure lb), which provides a 3<r upper 
limit to the time lag relative to V for 1.25 /im — 
25 /im of —50 days (phase — 0.16). The data for 
R Lep (Figure 2a) provide a 3<r upper limit to a 
possible 4.9 /im to V time lag of 20 days (phase — 
0.05), however, strong limits are not available for 
the other DIRBE wavelengths. 

4.4. Amplitude Variations with Wave- 
length 

Except for the phase lags noted above, in gen- 
eral the DIRBE light curves trace the optical 
light curves, but with smaller amplitudes of varia- 
tion. For example, the optical brightness of R Hor 
dropped 4.5 magnitudes between JD = 2448000 
and 2448160 (a factor of 63 in flux density), while 
the 1.25 micron flux density decreased by only a 
factor of 2.4 during this time (0.6 magnitudes). In 
general, the amplitude decreases with increasing 
wavelength. This is shown in Figure 4b, where 
the fitted amplitudes for the three stars in Table 
3 are plotted as a function of wavelength. 

As noted earlier, it was not possible to fit the 
light curves for most of our stars to an assumed si- 
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nusoidal variation, due to the incompleteness and 
peculiarities of the light curves. To investigate am- 
plitude variations in a larger sample of stars, for 
the stars with near-continuous high S /N time cov- 
erage for at least 10 weeks in both the optical and 
DIRBE data, in Figure 5a we plot the observed 
change in the V magnitude vs. the observed change 
in the 1.25 /im magnitude over the time period 
for which we have near-continuous data. We em- 
phasize that these magnitude changes generally do 
not represent the full range of variation for these 
stars: they only represent time ranges for which 
we have full sets of data. We have excluded time 
periods during which a maximum or minimum was 
occuring because of the possibility of phase lags. 
In Figure 5b, we plot Amag(V) against Amag(4.9 
jitm), using the same criteria. 

In Figures 5c — 5g, we plot Amag(1.25 /xm) vs. 
Amag(2.2 /xm), Amag(2.2 /xm) vs. Amag(3.5 /im), 
Amag(3.5 /im) vs. Amag(4.9 /im), Amag(4.9 /im) 
vs. Amag(12 /im), and Amag(12 /im) vs. Amag(25 /tm). 
The data plotted in these figures is for all stars in 
our sample, regardless of whether optical data are 
available, excluding only stars with less than 10 
weeks of high S /N data. The magnitude changes 
were measured between the observed maxima and 
minima of the light curves, ignoring time periods 
during which no data were available. As with the 
Figures 5a and 5b, these magnitude changes do 
not necessarily represent the full range of varia- 
tion for these stars; they are merely the magnitude 
changes observed in the DIRBE data. 

Figures 5a and 5b show that, although the V 
magnitudes of these stars may vary by several 
magnitudes during the plotted time intervals, the 
1.25 /im magnitudes typically only change by ~0.5 
magnitudes, and the 4.9 /xm typical vary by only 
^0.25 magnitudes. On average, for every mag- 
nitude of variation in V, the 1.25 /im brightness 
varies by 0.20 ± 0.01 magnitudes and the 4.9 
/xm brightness varies by 0.14 ± 0.01 magnitudes. 
There is a slight correlation of A (1.2 5 /im) with 
AV, and no observed correlation of A(4.9 /im) 
with AV. 

In all of the infrared-infrared plots (Figures 
5c-5g), with the exception of A(25 /im) vs. 
A (12 /im), clear correlations are visible. The 
slopes of the best-fit lines are all less than 
one, indicating decreasing amplitude, on average, 
with increasing wavelength. For each magnitude 



change at 1.25 /im, on average, the 2.2 /xm flux 
changes by only 0.69 ±0.01 magnitudes. The slope 
of the Amag(2.2 /im) vs. Amag(3.5 /xm) plot is 
0.83 ± 0.01, while the best-fit Amag(3.5 /im) vs. 
Amag(4.9 /im) line has a slope of 0.96 ± 0.01. 
The Amag(4.9 /xm) vs. Amag(12 /xm) best-fit line 
has a slope of 0.93 ± 0.04, and the Amag(12 /xm) 
vs. Amag(25 /im) slope is 0.92 ± 0.10. Thus 
the decrement of the amplitude with wavelength 
is smaller in the mid-infrared than in the near- 
infrared. 

In these plots, the slopes for oxygen-rich and 
carbon-rich stars are similar, except that oxygen- 
rich stars have somewhat larger Amag(2.2 /xm) /Amag(1.25 /xm) 
ratios on average (0.69 ± 0.01) than the carbon- 
rich stars (0.39 ± 0.04), and somewhat smaller 
Amag(4.9 /iin)/Amag(3.5 /xm) ratios (0.94 ± 
0.01 vs. 1.26 ± 0.05). There are no signifi- 
cant differences between the slopes of long pe- 
riod (> 400 days) and short period stars, or 
between early-type (M4 — M6.5) and late-type 
(M7 — M10) oxygen-rich stars. The slopes for 
stars with dust spectral types of SE1 — SE4 
are consistent with those for stars with types 
SE5 - SE7 (Table 2), with the exception of 
the Amag(12 /xm)/Amag(4.9 /im) ratio, which 
is higher in SE5 - SE7 stars (1.02 ± 0.02 vs. 0.71 
± 0.05). Spectral types SE5 - SE7 have nar- 
rower, more pronounced 9.7 /xm silicate emission 
features than types SE1 - SE4 (Sloan & Price 
1998). As noted earlier, a number of multi-epoch 
infrared spectroscopic studies show that, at least 
for some Miras, as the star increases in bright- 
ness, the strength of the silicate feature tends to 
increase (Little-Marenin & Stencel 1992; Hron & 
Aringer 1994; Little-Marenin et al. 1996; Creech- 
Eakman et al. 1997; Onaka et al. 1997; Monnier et 
al. 1998). Thus, over a pulsation period, for stars 
with strong silicate emission features (dust types 
SE5 — SE7), this effect may cause a larger change 
hi F12 fim/Fi.g than in other stars, enhancing 
the 12 /xm/4.9 /im amplitude ratios for these stars 
compared to stars with dust types SE1 — SE4. 

Amplitude decrements with wavelength are also 
visible in some of the spectral energy distribution 
plots (Figure 3). In particular, note the variations 
in the spectral energy distributions of o Cet and 
IK Tau with time (Figures 3c and 3d). At a given 
wavelength, the separation between data points 
obtained at different times tends to get smaller as 
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the wavelength increases. 

4.5. General Remarks about the Optical- 
Infrared Spectra 

Figure 3 shows that the variation in the spectral 
energy distribution of a given star over a pulsation 
period is much less than star-to-star differences 
in spectral shape. The majority of the stars in 
our sample have spectral energy distributions that 
peak at 2.2 /um, but several are brightest at 1.25 
/zm (for example, X Oph and R Aql; Figure 3h), 
and a number of stars peak at 3.5 /zm (for exam- 
ple, R For and V Cyg; Figure 3a). Three stars, LP 
And (Figure 3b), V384 Per (Figure 3a), and WX 
Psc (Figure 3c), have very red spectral energy dis- 
tributions, peaking at 4.9 /tzm. Figure 3 shows that 
there is a clear statistical difference in the infrared- 
optical spectral energy distributions of carbon and 
oxygen- rich Miras, with carbon stars tending to 
have redder spectral energy distributions. Six of 
the 8 carbon stars peak at 3.5 /jm or 4.9 /zm, but 
only two of the 30 oxygen stars. 

We note that WX Psc, the only oxygen-rich star 
in our sample to peak at 4.9 (im (Figure 3c), is 
also the only star in our sample with an IRAS 
dust spectral type of 'silicate self-absorbed emis- 
sion' (SB) (Table 2; see Sloan & Price (1998)). 
The rest of the stars have silicate emission spectra 
(dust class SE). With the exception of WX Psc, 
within the oxygen-rich class there is no strong cor- 
relation of the IRAS dust spectral type (Table 2) 
with the peak of the infrared-optical spectral en- 
ergy distribution as seen in the DIRBE data (Fig- 
ure 3). There is a weak correlation between opti- 
cal spectral type (Table 2) and the wavelength of 
the peak of the optical-infrared spectrum (Figure 
3), in that the two oxygen-rich stars with infrared 
spectral energy distributions that peak at 3.5 /zm 
(IW Hya; Figure 3f) and 4.9 /im (WX Psc; Figure 
3c) have very late optical spectral types (M9 and 
M8, respectively). 

For the carbon-rich stars, all but two have IRAS 
dust spectral types of SiC, i.e., with the 11.2 /zm 
SiC feature (Table 1, see Sloan et al. (1998)). One 
star (S Cep) has an IRAS spectral type of SiC+ 
(the 11.2 /zm feature plus a secondary feature at 
8.5 — 9 /an). The last carbon star in our sam- 
ple, LP And, has an IRAS classification of 'red', 
meaning that, in addition to the 11.2 /zm feature, 
in the IRAS 8 — 23 /zm spectrum the star shows a 



strong infrared continuum increasing with increas- 
ing wavelengths. LP And is one of the two carbon 
stars in our sample with a very red DIRBE spec- 
trum, peaking at 4.9 /jm (Figure 3b). LP And is 
also the carbon star with the latest optical spectral 
type in our sample, C8. For the rest of the car- 
bon stars, excluding LP And, there is no strong 
correlation of optical spectral type (Table 1) with 
DIRBE spectral peak (Figure 3). 

5. Discussion 

Our data show that the amplitude of varia- 
tion for our sample stars decreases with increasing 
wavelength, from the optical out to 25 /tni. The 
amount of decrease with wavelength is consistent 
with the ground-based results of Le Bertre (1992, 
1993). If, at a particular time during the pulsation 
cycle, one approximates the ratio of the fluxes in 
two adjacent wavelength bands by an attenuated 
blackbody spectrum, the observed amplitude de- 
crease with increasing wavelength means that the 
temperature of this blackbody increases as the star 
becomes brighter. 

The observed amplitude decrement with wave- 
length of a particular carbon-rich Mira, R For, was 
successfully modeled by Le Bertre (1988). In this 
model, the dust condensation temperature is fixed, 
which causes the radius of condensation to move 
in and out in phase with the luminosity of the 
star. The stellar spectrum was approximated by a 
blackbody, and a A -13 dust opacity law was used. 
At a given radius, the dust temperature is high- 
est at highest stellar luminosity. Averaged over all 
the dust in the shell, the effective dust tempera- 
ture is highest at maximum, consistent with the 
observations. 

More general, and more complex, models for 
carbon-rich long period variables were developed 
by Winters et al. (1994b, 2000), using time- 
dependent hydrodynamics and a detailed treat- 
ment of the formation, growth, and evaporation 
of carbon grains. In these models, the pulsation 
of the underlying star is simulated by a sinusoidal 
variation of the velocity in the lowest layer of the 
atmosphere, and isothermal shocks and complete 
dust-grain momentum coupling are assumed. New 
layers of dust are formed periodically in the cir- 
cumstcllar envelope, and radiation pressure pushes 
this dust, and its swept-up gas, outwards. Dust 



8 



formation and growth are triggered by compres- 
sion of the gas due to the pulsation, and grain nu- 
cleation, growth, and evaporation are calculated 
according to the relations in Gail, Keller & Sedl- 
mayr (1984), Gail & Sedlmayr (1988), and Gauger 
ct al. (1990). The frequency-dependent radiative 
transfer equation in spherical geometry is solved as 
in Winters, Dominik, & Sedlmayr (1994a), using 
gray absorption and thermal emission of the gas 
and frequency-dependent absorption and thermal 
emission of the dust. In the calculation of the dust 
opacity, the particles are assumed to be spherical 
and small compared to the relevant wavelength of 
the radiation field. 

Using this model, Winters et al. (1994b) calcu- 
lated light curves for a 1 M Q , 10 4 L Q star with an 
assumed stellar temperature of 2600K, a period of 
650 days, and a velocity amplitude of 2 km s _1 . 
They find amplitudes of 8 magnitudes at 0.55 /im, 
2.8 magnitudes at 0.9 /im, 1 magnitude at 1.65 
^m, 0.7 magnitudes at 2.2 /im, 0.35 magnitudes 
at 4.8 /on, and 0.2 magnitudes at 25 /im. These 
amplitude ratios are within a factor of two of our 
average values for carbon stars, with the exception 
of the 25 /im amplitude, which is low by a factor 
of four. 

By comparing the DIRBE data with light 
curves from the AAVSO, we found evidence for 
time lags in the near-infrared maxima (1.25 — 3.5 
/im) relative to the optical for six stars, with phase 
shifts between ~0.05 — 0.13. This confirms that 
the stellar luminosity maximum occurs after visual 
maximum, and the V light curve does not trace 
the bolomctric luminosity. We also found possible 
lags between the near-infrared and mid-infrared 
maxima (4.9 — 12 /im) for three stars, with the 
mid-infrared flux peaking before the near-infrared 
but after visual maximum. The mid-infrared to 
near-infrared relative phase shifts are between ~ 
0.05 - 0.12. 

At present, there are no theoretical models of 
Miras that predict both of these lags. The well- 
known 1.04 /im phase lag in Miras relative to V 
(Lockwood & Wing 1971; Barnes 1973) has been 
successfully modeled by Alvarez & Plez (1998) 
as a consequence of the propagation of perturba- 
tions through the atmosphere. These perturba- 
tions cause time lags in the strength of the vana- 
dium oxide absorption relative to titanium ox- 
ide, because these species are formed at different 



depths in the atmosphere. The V magnitudes of 
M stars are strongly attenuated by TiO absorption 
(Smak & Wing 1979a). 

Another factor that may contribute to phase 
lags is variations in the optical depth of the cir- 
cumstellar shell, which may periodically attenu- 
ate the visual light more than the near- and mid- 
infrared. Extinction-induced infrared-optical off- 
sets in the maxima are seen in the models of Win- 
ters et al. (1994b), however, in these models, op- 
tical maximum occurs after infrared maximum, 
while in our observations, the infrared maximum 
lags that in the optical. We note that the Winters 
et al. (1994b) models do not include the effect of 
molecular line opacity, while the models of Alvarez 
& Plez (1998) do not include dust absorption, so 
neither set of models do a complete analysis of 
phase lags. It is possible that, for some stars, dust 
extinction can be neglected, for example, Lobel et 
al. (2000) model variations in the spectral energy 
distribution of o Cct using the DUSTY radiative 
transfer code (Ivezic, Nenkova, & Elitzur 1996), 
and find that the circumnuclear shell of this star 
is optically thin. They conclude that the observed 
visual-to-near-infrared phase lag of Mira cannot 
be due to modulation of the optical by circumstel- 
lar extinction, but instead must originate in the 
stellar atmosphere itself. 

In Section 3.2, we noted that some of our light 
curves are asymmetric, with secondary maxima or 
plateaus in the rising portion of the light curve. 
These two types of light curves were previously 
noted by Le Bertre (1992). Similar features ap- 
pear in the models of Winters et al. (1994b, 1995), 
where they are due to the formation of a new 
dust layer and subsequent excess infrared dust 
emission. Alternative explanations for such inflec- 
tion points are dissipation of shock wave energy 
in the inner part of the atmosphere (Fcuchtingcr 
et al. 1993) or formation of secondary shocks in 
the inner-most dust-free region (Bowen et al. 1990; 
Willson et al. 1994). The fact that matching opti- 
cal excesses are present in the AAVSO data sup- 
ports the hypothesis that these features are due 
to shocks rather than dust, since the dust models 
predict that extinction by the newly-formed dust 
layer will cause concurrent depressions (rather 
than enhancements) of the optical light curves. 
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6. Conclusions 

Using data from the DIRBE instrument on 
COBE, we have obtained 1.25 /im — 25 fim light 
curves for 38 infrared-bright Mira variables and 
compared with AAVSO optical data. We have 
found decreasing amplitudes of variation with in- 
creasing wavelength, with magnitudes consistent 
with recent theoretical models of circumstcllar 
dust shells around Mira stars. For the oxygen- 
rich Miras R Car and T Cep and the carbon star 
S Cep, we see secondary maxima in the rising part 
of the DIRBE light curves. These features are also 
present in the optical data from the AAVSO, sug- 
gesting that they are due to shocks rather than to 
dust formation. 

We have detected near-infrared-to-optical time 
delays of the maxima of six stars with phases ~ 
0.05 to 0.12. We also find possible offsets of the 
mid-infrared maxima relative to those in the near- 
infrared for three stars, in that the mid-infrared 
peaks occurs after optical maximum but before 
that in the near-infrared. These lags may be 
caused by variations in the amount of molecular 
absorption over a pulsation period, in combina- 
tion with variations in dust opacity. More detailed 
modeling is needed to disentangle these two ef- 
fects. 

We find no evidence for large magnitude (>0.2 
magnitudes) short-term (hours to days) variations 
in the near- or mid-infrared in our stars. With 
typical time resolution of 9 hours, however, the 
COBE data are not particularly sensitive to bursts 
of a few hours in duration. 
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used, and acknowledge with thanks, data from the 
AAVSO International Database, based on obser- 
vations submitted to the AAVSO by variable star 



observers worldwide. This research has made use 
of the SIMBAD database, operated at CDS, Stras- 
bourg, France. This research was funded in part 
by NSF POWRE grant AST-0073853. 

REFERENCES 

Alvarez, R., & Plcz, B. 1998, A&A, 330, 1109 

Anandarao, B. G., Pottasch, S. R., & Vaidya, D. 
B. 1993, A&A, 273, 570 

Barnes, T. G. 1973, ApJS, 25, 369 

Barthes, D., Chcncvez, J., & Mattei, J. A. 1996, 
AJ, 111, 2391 

Bastian, U. & Roser, S. 1991, Sterne und Wel- 
traum, 30, 592 

Bedijn, P. J. 1987, A&A, 186, 136 

Boggess, N. W., et al. 1992, ApJ, 397, 420 

Bowen, G. H. 1990, in The Numerical Modelling 
of Non-linear Stellar Pulsations, Buchler, J. R. 
(ed), Kluwer Academic Publishers, p. 155 

Catchpole, R. M., Robertson, B. S. C, Lloyd- 
Evans, T. H. H., Feast, M. W., Glass, I S., & 
Carter, B. S. 1979, SAAO Circulars, 1, 61 

Claussen, M. J., Klcinmann, S. G., Joyce, R. R., 
& Jura, M. 1987, ApJS, 65, 385 

Conrow, T., et al. 1993, B.A.A.S., 183, 03.03 

Creech-Eakman, M. J., Stencel, R. E., Williams, 
W. J., & Klcbe, D. I. 1997, ApJ, 477, 825 

De Laverny, P., Mcnnessier, M. O., Mignard, F., 
& Mattei, J. A. 1998, A&A, 330, 169 

Feuchtinger, M. U., Dorfi, E. A., & Hofncr, S. 
1993, A&A, 273, 513 

Gail, H.-P., Keller, R., & Sedlmayr, E. 1984, A&A, 
235, 245 

Gail, H.-P., & Sedlmayr, E. 1988, A&A, 206, 153 

Gaugcr, A., Gail, H.-P., & Sedlmayr, E. 1990, 
A&A, 235, 345 

Guenter, H., & Henson, G. 2001, in preparation 

Harvey, P. M., Bechis, K. P., Wilson, W. J., & 
Ball, J. A. 1974, ApJS, 27, 331 



10 



Hauser, M. G., Kelsall, T., Leisawitz, D., & Wei- 
land, J. 1998, 'COBE Diffuse Infrared Back- 
ground Experiment (DIRBE) Explanatory Sup- 
plement', V. 2.3. 

Hron, J. & Aringcr, B. 1994, A&A, 281, 139 

Ivezic, Z., Nenkova, M., & Elitzur, M. 1996, User 
Manual for DUSTY, available on the Internet: 
http://www.pa.uky.edu/moshc/dusty/ 

Ivezic, Z, & Knapp, G. R. 2000, Proceedings of 
IAU Symposium 191: AGB Stars, A.S.P. Con- 
ference Series, V. 203, p. 124 

IRAS Point Source Catalog (1988), Version 2, 
Joint IRAS Science Working Group (Washing- 
ton: GPO). 

Jarrett, T. 1992, Ph.D. Thesis, University of Mas- 
sachusetts. 

Jura, M. 1986, ApJ, 303, 327 

Knapp, G. R., Phillips, T. G., Leighton, R. B., 
Lo, K. Y., Wannicr, P. G., Woottcn, H. A., & 
Huggins, P. J. 1982, ApJ, 252, 616 

Knapp, G. R., & Morris, M. 1985, ApJ, 292, 640 

Koornncef, J. 1983, A&A, 128, 84. 

Kholopov et al. 1985-1988, General Catalog of 
Variable Stars. 

Le Bertre, T. 1988, A&A, 190, 79 

Le Bertre, T. 1992, A&AS, 94, 377 

Le Bertre, T. 1993, A&AS, 97, 729 

Le Bertre, T. 1993, A&A, 324, 1059 

Le Bertre, T., & Winters, J. M. 1998, A&A, 334, 
173 

Le Sidancr, P., & Bertre, T. 1993, A&A, 278, 167 

Le Sidaner, P., & Bertre, T. 1996, A&A, 314, 896 

Littlc-Marcnin, I. R. & Stencel, R. E. 1992, 
PASPC, 26, 591 

Littlc-Marcnin, I. R., Stencel, R. E., & Staley, S. 
B. 1996, ApJ, 467, 806 

Lockwood, G. W. & Wing, R. F. 1971, ApJ, 169, 
63 



Lobel, A., Bagnulo, S., Doyle, J. G., & Power, C. 
2000, MNRAS, 317, 391 

Lobel, A., Doyle, J. G., & Bagnulo, S. 1999, A&A, 
343, 466 

Maffei, P., & Tosti, G. 1995, AJ, 109, 2652 

Maran, S. P., ct al. 1977, Infrared Phys., 17, 565 

Mattei, J. A. 2000, Observations from the AAVSO 
International Database, private communica- 
tion. 

Monnier, J. D., Geballe, T. R., & Danchi, W. C. 
1998, ApJ, 502, 833 

Moshir, M., et al. 1992, Explanatory Supplement 
to the IRS A Faint Source Survey, Version 2, 
JPL D-10015 8/92 (Pasadena: JPL) 

Murakami, H., et al. 1996, P.A.S.J., 48, L41 

Olofsson, H. 1999, Proceedings of the Science with 
the Atacama Large Millimeter Array, 45. 

Onaka, T., de Jong, T., & Willems, F. J. 1989a, 
A&AS, 81, 261 

Onaka, T., dc Jong, T., & Willems, F. J. 1989b, 
A&A, 218, 169 

Onaka, T. et al. 1997, Ap&SS, 255,311 

Petit, E., & Nicholson, S. B. 1933, ApJ, 78, 320 

Sedlmayr, E. 1994, in Molecules in the Stellar 
Environment, IAU Colloquium 146, (Springer- 
Verlag: New York), 163 

Skrutskie, VI., et al. 1997, in The Impact of Large- 
Scale Near-IR Sky Surveys, ed. F. Garzon et al. 
(Dordrecht: Kluwer), 187 

Sloan, G. C, & Price, S. D. 1995, ApJ, 451, 758 

Sloan, G. C, & Price, S. D. 1998, ApJS, 119, 141 

Sloan, G. C, Little-Marenin, I. R., & Price, S. D. 
1998, AJ, 115, 809 

Smak, J., & Wing, R. F. 1979, Acta. Astron., 29, 
187 

Smak, J., & Wing, R. F. 1979, Acta. Astron., 29, 
199 



11 



Smith, B. J. 1995, J Band Catalog 
Information, Internal IPAC Memo, 

http://www.etsu.edu/physics/bsmith/irtsdat/JBAND.CAT.VER5.INFO 

Spear, G. C, Owen, B., & Madore, B. 1993, 
PASPC, 43, 35 

Wheaton, W. A., et al. 1997, Pointing Reconstruc- 
tion Processing for IRTS, Internal IPAC Memo, 
http://spider.ipac.caltech.edu/staff/waw/irts/irts.html 

Willson, L. A., Struck-Marcell, O, & Bowen, G. 
H. in Cosmic Winds and the Heliosphcre, eds. 
J. R. Jokipii, C. P. Sonett, and M. S. Giampapa 

Winters, J. M., Dominik, C, & Sedlmayr, E. 
1994a, A&A, 288, 255 

Winters, J. M., Fleischer, A. J., Gauger, A., & 
Sedlmayr, E. 1994, A&A, 290, 623 

Winters, J. M., Fleischer, A. J., Gauger, A., & 
Sedlmayr, E. 1995, A&A, 302, 483 

Winters, J. M., Le Bertre, T., Jeong, K. S., 
Helling, C, & Sedlmayr, E. 2000, A&A, 361, 
641 

Wood, P. R. 1990, in Confrontation between Stel- 
lar Pulsations and Theory, eds. C. Cacciari & 
G. Clementini, ASP Conf. Series 11, 355 

Wood, P. R. 1997, in Planetary Nebulae, IAU 
Symp. 180, 297 

Woolf, N. J., & Ncy, E. P. 1969, ApJ, 155, L181 



This 2-column preprint was prepared with the AAS IATjrjX 
macros v5.0. 



12 



Fig. 1.— The COBE DIRBE 1.25 - 25 fim light 
curves from the DIRBE Weekly Maps for six sam- 
ple stars with relatively complete light curves: S 
Cep, R Hor, RW Vol, R Car, T Cep, and R Cas. 
The carbon star S Cep is plotted first, followed by 
the five oxygen-rich stars in R.A. order. When 
available, optical light curves from the AAVSO 
are also plotted. The 12 fim/ 2. 2 /im flux ratio 
as a function of time is also provided. For RW 
Vel, which does not have optical data available, 
F4.9 Aim /F2.2 A «n is plotted instead of the optical 
light curve. Note that the infrared light curves 
are plotted in flux density units, while the optical 
data are given in magnitudes. 

Fig. 2.— The 4.9 ^m DIRBE light curves for the 
remaining 31 stars in our sample, from the DIRBE 
Weekly Maps. The carbon stars (Table 1) are plot- 
ted first, in R.A. order, followed by the oxygen-rich 
stars (Table 2). These are compared with optical 
light curves when possible. 

Fig. 3. — The optical to infrared spectra for se- 
lected stars, at one; to three different times during 
the COBE mission. The stars are plotted in R.A. 
order, with the carbon stars first, as in Tables 1 
and 2. The symbols indicate the location in the 
light curve as follows: filled triangles: times near 
minimum; open triangles: inflection points; open 
circles: near infrared maximum; filled circles: op- 
tical maximum (when clearly offset from infrared 
maximum); open squares: times when the light 
curve is rising; crosses: times when the light curve 
is falling. A hexagon is used when it is unclear 
what part of the light curve the data are located. 
All data points within 1 day of the given time were 
plotted. Note that the shape of the spectra vary 
with time, with the shorter wavelength fluxes in- 
creasing a larger amount with increasing bright- 
ness than longer wavelength fluxes. Also note that 
the spectra of carbon stars tends to peak at longer 
wavelengths than those of oxygen-rich stars. 



Table 3. The symbols are the same as those in 
Figure 4a. 

Fig. 5. — a) Top left: Amag(1.25 fim) vs. 
Amag(V) for the stars with at least 10 continu- 
ous weeks of optical and infrared data available. 
The uncertainties in V were calculated assum- 
ing ^0.5 magnitudes uncertainties in the AAVSO 
data. Note that the plotted magnitude changes do 
not always represent the full amplitude of the star; 
they are merely the observed change in magnitude 
over the available continuous data, b) Top right: 
Amag(4.9 fim) vs. Amag(V) for the stars with at 
least 10 continuous weeks of optical and infrared 
data available, c) Second row, left: Amag(1.25 
fim) vs. Amag(2.2 fim) for the stars with at least 
10 weeks of high S/N data at both wavelengths. 
The magnitude changes were measured between 
the observed maxima and minima of the light 
curves, ignoring time periods during which no data 
were available. Note that the plotted magnitude 
changes do not always represent the full amplitude 
of the star; they are merely the observed change 
in magnitude over the available DIRBE data (ig- 
noring gaps in the light curves), d) Second row, 
right: Amag(2.2 /im) vs. Amag(3.5 fim) for the 
stars with at least 10 weeks of high S/N data at 
both wavelengths, e) Third row, left: Amag(3.5 
fim) vs. Amag(4.9 fim) for the stars with at least 
10 weeks of high S/N data at both wavelengths, 
f) Third row, right: Amag(4.9 /jm) vs. Amag(12 
fim) for the stars with at least 10 weeks of high 
S/N data at both wavelengths, g) Fourth row, 
left: Amag(12 fim) vs. Amag(25 fim) for the stars 
with at least 10 weeks of high S/N data at both 
wavelengths. 



Fig. 4. — a) Phase lag of maximum relative to 
the maximum at V vs. wavelength, for the three 
stars with sine fits to their light curves (Table 3): 
R Hor (filled triangles), R Cas (open circles), and 
RW Vel (crosses). No appropriate optical data 
were available for RW Vel, so we assumed a 1.25 
fim lag relative to V of 0.1 phase, b) Amplitude 
of variation vs. wavelength, for the three stars in 
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Table 1 

Carbon-Rich Mira Variable Stars in Sample 



Name 


Period 


[12] 


Spectral 


Dust 


Optical 




(Days) 


(mag) 


Type 


Spectrum 


Data? 


R For 


388.73 


-2.38 


C4,3e(Ne) 


SiC 


Yes 


V384 Per 


535.0 


-3.19 


C(N) 


SiC 




R Lep 


427.07 


-2.82 


C7,6e(N6e) 


SiC 


Yes 


VCyg 


421.27 


-3.43 


C5,3e-C7,4e(Npe) 


SiC 


Yes 


RV Aqr 


453.50 


-2.59 


C6-7,2-4(Ne) 


SiC 




V1426 Cyg 


470.0 


-2.40 


C7,2e(N) 


SiC 




S Cep 


486.84 


-2.83 


C7.4e,(N8e) 


SiC+: 


Yes 


LP And 




-3.83 


C8,3.5e 


Red 
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Table 2 

Oxygen-Rich Mira Variable Stars in Sample 



Name 


Period 


[12] 


Spectral 


Dust 


Optical 




(Days) 


(mag) 


Type 


Spectrum 


Data? 


KU And 


750.00 


-3.06 


M10 I-III 


SE5 




T Cas 


444.83 


-2.95 


M6-9.0e 


SE1 


Yes 


WX Psc 


660.00 


-4.03 


M8 


SB 




o Cet 


331.96 


-5.59 


M5-9e 


SE8 


Yes 


R Hor 


407.60 


-3.53 


M5-8eII-III 


SE5t 


Yes 


IK Tau 


470.00 


-5.54 


M6-10e 


SE7 


Yes 


TX Cam 


557.40 


-4.41 


M8-10 


SE6 


Yes 


R Aur 


457.51 


-3.03 


M6.5-9.5e 


SE2 


Yes 


U Ori 


368.30 


-3.46 


M6-9.5e 


SE6t 


Yes 


GX Mon 


527.00 


-3.32 


M9 


SE6 




R Cnc 


361.60 


-2.54 


M6-9e 


SE2 


Yes 


XYVV Vei 






M7 TTlYTTV 

1V1 / 111 111 1 L. 






R Car 


308.71 


-3.06 


M4-8e 


SE1 


Yes 


R LMi 


372.19 


-2.94 


M6.5-9.0e(Tc:) 


SE4t 


Yes 


IW Hya 


650.00 


-3.32 


M9 


SE5 




R Leo 


309.95 


-4.71 


M6-8 III-9.5e 


SE2 


Yes 


AQ Cen 


387.50 


-2.32 


Me 


SE7 




R Cen 


546.20 


-3.09 


M4-8IIe 


SE2 


Yes 


WX Ser 


425.10 


-2.30 


M8e 


SE6 




U Her 


406.10 


-3.12 


M6.5-9.5e 


SE4 


Yes 


Villi Oph 




-3.51 


M4III-9 


SE5 




X Oph 


328.85 


-2.90 


M5-9e 


SE1 


Yes 


V3953 Sgr 




-3.39 


M9 


SE6 




R Aql 


284.20 


-2.88 


M5-9e 


SE5 


Yes 


V3880 Sgr 


510.00 


-2.53 


M8: 


SE5 




V342 Sgr 


372.00 


-2.63 


M9 


SE6 




RR Aql 


394.78 


-2.67 


M6e-9 


SE7 


Yes 


T Cep 


388.14 


-3.56 


M5.5-8.8e 


SE1 


Yes 


R Aqr 


386.96 


-4.37 


M5-8.5e+pcC 


SE7 


Yes 


R Cas 


430.46 


-4.19 


M6-10e 


SE5t 


Yes 
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Table 3 

Parameters of Best-Fit Sine Curves 



Star 


Wavelength 


Date of Maximum 


Amplitude 




W 


(JD) 


(Magnitudes) 


R Hor 


0.55 


2447934 


7.66 




1.25 


2447987 


1.14 




2.2 


2447999 


0.83 




3.5 


2447993 


0.56 




4.9 


2447959 


0.45 




12 


2447939 


0.59 




25 


2447960 


0.64 


R Cas 


0.55 


2447957 


5.53 




1.25 


2447994 


1.42 




2.2 


2448003 


1.03 




3.5 


2447995 


0.84 




4.9 


2447982 


0.77 




12 


2447978 


0.70 




25 


2447980 


0.68 


RW Vel 


1.25 


2448013 


0.80 




2.2 


2448032 


0.59 




3.5 


2448025 


0.41 




4.9 


2447983 


0.29 
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